# Metview Macro

temperature_forecast = read("temperature_forecast.grib")
temperature_analysis = read("temperature_analysis.grib")

print(type(temperature_forecast))
print(count(temperature_forecast))

fc_an_diff = temperature_forecast - temperature_analysis

neg_shade = mcont(
    legend                         : "on",
    contour                        : "off",
    contour_highlight              : "off",
    contour_level_selection_type   : "level_list",
    contour_level_list             : [-30,-10,-5,-3,-1.5,-0.5],
    contour_label                  : "off",
    contour_shade                  : "on",
    contour_shade_method           : "area_fill",
    contour_shade_max_level_colour : "cyan",
    contour_shade_min_level_colour : "blue",
    contour_shade_colour_direction : "clockwise"
    )

pos_shade = mcont(
    legend                         : "on",
    contour                        : "off",
    contour_highlight              : "off",
    contour_level_selection_type   : "level_list",
    contour_level_list             : [0.5,1.5,3,5,10,30],
    contour_label                  : "off",
    contour_shade                  : "on",
    contour_shade_method           : "area_fill",
    contour_shade_max_level_colour : "red",
    contour_shade_min_level_colour : "yellow",
    contour_shade_colour_direction : "clockwise"
    )

plot(max(fc_an_diff), pos_shade, neg_shade)
#plot(min(fc_an_diff), pos_shade, neg_shade)
#plot(mean(fc_an_diff), pos_shade, neg_shade)
